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In dense hot star winds, the infrared and radio continua are dominated by free-free opacity and recombination 
emission line spectra. In the case of a spherically symmetric outflow that is isothermal and expanding at constant 
radial speed, the radiative transfer for the continuum emission from a dense wind is analytic. Even the emission 
profile shape for a recombination line can be derived. Key to these derivations is that the opacity scales with only 
the square of the density. These results are well-known. Here an extension of the derivation is developed that also 
allows for line blends and the inclusion of an additional power-law dependence beyond just the density dependence. 
The additional power-law is promoted as a representation of a radius dependent clumping factor. It is shown that 
differences in the line widths and equivalent widths of the emission lines depend on the steepness of the clumping 
power-law. Assuming relative level populations in LTE in the upper levels of He II, an illustrative application of the 
model to SpitzerflRS spectral data of the carbon-rich star WR 90 is given. 
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1 Introduction 

The winds of massive stars are exceedingly dense, espe- 
cially in evolved stars such as OB supergiants, Luminous 
Blue Variables (LBVs), and Wolf-Rayet (WR) stars (e.g., 
Cassinelli 1979; Lamers & Cassinelli 1999; Kudritzki & 
Puis 2000). The winds are dense and highly ionized, to the 
point where free-free opacity can become optically thick 
in the wind not only at radio wavelengths but also in the 
infrared (IR). Consequently, the IR band offers opportuni- 
ties to study the wind density distribution through contin- 
uum and line emission. Wright & Barlow (1975) showed 
that at radio wavelengths the specific flux of emission for 
a spherical, constant expansion, isothermal, and optically 
thick wind will produce a continuum power-law having a 
slope of f v oc ^ 6 , where /„ is the specific flux (e.g., in 
Janskys). Knowing the distance to the star, the specific lu- 
minosity is derivable and can be related to the wind mass- 
loss rate M. Observations of appropriate sources at radio 
frequencies have been important for deriving M values 
(e.g., Abbott et al. 1980; Leitherer, Chapman, & Koribal- 
ski 1995). 

Understanding mass loss during the various evolution- 
ary stages of massive stars is of critical importance, as for 
example in understanding observed ratios of blue and red 
supergiant stars or the relative numbers of WR subtypes, 
the nitrogen, carbon, and oxygen rich WR stars (Maeder 
& Meynet 2000; Meynet & Maeder 2003). Mass-loss rates 
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are also important for interpreting the afterglows of gamma- 
ray bursts, of which some are associated with supernova 
explosions (Woosley & Bloom 2006). 

Abbott, Bieging, & Churchwell (1981) showed that 
radio emission used to infer mass-loss rates can also de- 
pend on "clumping" if stellar winds are not laminar or 
"smooth". The reason is that the free-free opacity scales 
with the square of the wind density. As a result, more 
or less clumping biases the mass-loss rate determinations 
if clumping corrections are not taken into account, with 
ramifications for our understanding of massive star evo- 
lution. And not only will clumping affect mass-loss rates, 
it can even have feedback for the line-driving efficiency 
of the winds in some cases (Brown et al. 2004; Oskinova, 
Hamann, & Feldmeier 2007). 

It has become abundantly clear that stellar winds are 
indeed clumped, and that the clumping is not negligible 
in relation to interpreting M values. Already, researchers 
had good reason to expect clumping to be important, since 
the hot star winds are understood to be line-driven (Cas- 
tor, Abbott, & Klein 1975; Friend & Abbott 1986; Paul- 
drach, Puis, & Kudritzki 1986), a mechanism that is sub- 
ject to instability leading to the formation of wind struc- 
ture (Lucy & Solomon 1970; Owocki, Castor, & Rybicki 
1988; Feldmeier, Puis, & Pauldrach 1997). Evidence for 
clumping in WR star winds and now even O star winds 
has been mounting for many years. In the case of WR 
stars, a convincing argument was made by Hillier (1991) 
on the basis of the emission of recombination line cores 
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in relation to the strength of the electron scattering wings 
of the lines. The former is a density square process and 
subject to influence by clumping, but the latter is linear 
in density. Smooth wind models were found inadequate 
in matching the line profile shapes, but the inclusion of 
clumping allowed for superior matches between models 
and data. 

Information about the clumping in the less dense winds 
of O stars has come only more recently (Bouret et al. 
2003; Evans et al. 2004; Bouret, Lanz, & Hillier 2005). 
Indeed, quite startling results have been reported by Fuller- 
ton, Massa, & Prinja (2006) who suggest that on the basis 
of FUV doublet lines of P V that early O star mass-loss 
rates might need downward revisions by factors of 10 to 
perhaps 100. This result along with the general recogni- 
tion of the prevalence of clumping in massive star winds 
led to an entire meeting devoted to the topic (Hamann, 
Feldmeier, & Oskinova 2007). 

Many researchers have shifted attention from the ques- 
tion of whether the wind clumping exists to its origin 
and evolution (i.e., how clumping initiates, how it evolves 
through the wind, and how its nature changes with stellar 
evolutionary phase). A better understanding of clumping 
requires fresh considerations of what diagnostics may be 
used to constrain clump properties, such as size and mass 
distribution and dynamics. For example, Nugis, Crowther, 
& Willis (1998) considered how a gradient in the clump- 
ing factor with radius in the wind modifies continuum 
slopes for the wind emission. The main point of their 
study is that a constant clumping factor affects the amount 
of emission, but it does not change the observed power- 
law slope of the continuum. Only a clumping factor that 
is radius dependent will influence the slope. 

Using the results of hydrodynamic simulations, Runacn 
& Owocki (2002, 2005) have considered how a structured 
flow evolves in radius for 1 -dimensional wind models. 
Necessarily, such models predict clumping in the form 
of spherical shells. This has long known not to be an 
accurate description of hot star winds, for example ow- 
ing to the relatively low variability observed in the X- 
ray emissions from massive stars (Cassinelli & Swank 
1983; Berghofer 1997). Indeed, there have been attempts 
to model the low X-ray variability in terms of statistically 
spherical but 3-dimensionally structured wind flow (Oski- 
nova et al. 2001). Although impressive progress has been 
made in extending the 1 -dimensional numerical simula- 
tions to an axisymmetric geometry (Dessart & Owocki 
2002, 2003, 2005), the problem of fully 3-dimensional 
and time-dependent radiative hydrodynamic models re- 
mains a challenge. 

With so much emphasis being given to measuring wind 
clumping factors and understanding structured wind flow, 
useful insight and tools can be gained from considering 
some simplified cases. As was noted, Wright & Barlow 



(1975) were able to derive a power-law continuum slope 
of index /„ oc at radio wavelengths for a spherical, 
isothermal, and inverse square law density. Still assuming 
a smooth wind flow, Cassinelli & Hartmann (1977) ex- 
tended the approach to allow for power-law density and 
temperature distributions. They derived analytic expres- 
sions for the continuum slope. And indeed, IR spectra of 
WR winds are known to be steeper than the canonical i/ - 6 
rule (e.g., Morris et al. 1993), and so approaches like that 
of Cassinelli & Hartmann are relevant. 

Emission line profile shapes have also been consid- 
ered. Hillier, Jones, & Hyland (1983) derived an analytic 
solution for the radiative transfer of the line and contin- 
uum emission for density-squared opacities. For the con- 
tinuum part, the standard Wright & Barlow (1975) solu- 
tion is recovered in terms of the dependence on the wind 
parameters and the power-law scaling in frequency. The 
inclusion of continuous opacity and line optical depth in 
the Sobolev approximation yields a line profile shape that 
is centrally peaked and gently rounded in appearance (see 
Ignace, Quigley, & Cassinelli 2003). 

Here the Hillier et al. (1983) solution is extended in 
two significant ways. First, their derivation prescribes a 
r~ 2 dependence for the wind density, but an additional 
power-law factor in radius can also be included in the 
opacity. Such a power law may represent the effects of 
a radius-dependent clumping factor. Second, it turns out 
that the radiative transfer can allow for an arbitrary num- 
ber of line blends, which is important for IR spectra dom- 
inated by emission lines. The result of the new formula- 
tion is an analytic solution to a fairly complex scenario 
of line blends and power-law opacity but in the restric- 
tive limit of spherical symmetry and constant expansion. 
Although direct applications may be limited, the solution 
provides complex cases against which sophisticated ra- 
diative transfer routines can be benchmarked. 

The following section describes the approach to the 
analytic solution that extends the results of Hillier et al. 
(1983), following the notation of Ignace et al. (2003). Sec- 
tion [3] provides an illustrative application to high resolu- 
tion spectra of the carbon-rich WC star, WR 90. A sum- 
mary of results is given in sectionS] 

2 Radiative Transfer Solution 
2.1 The Physical Ingredients 

The scenario being addressed assumes spherical symme- 
try for a wind expanding at constant radial speed. For a 
fixed mass-loss rate, this means that the density is an in- 
verse square law, hence 

4.7T Vac r 2 
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Only continuum and line processes that scale with the 
square of the density will be considered. This is a rea- 
sonable approximation for dense, ionized winds of hot 
stars where free-free continuum opacity and recombina- 
tion line opacity are important. The wind will be assumed 
isothermal. For the line transfer, the Sobolev approxima- 
tion will be adopted. 

Keys to the continuum radiative transfer are the ex- 
pressions for the continuum and line optical depths. For 
free-free opacity, the optical depth from the observer to a 
point in the wind is 



/>oo 

t v ,s = R* Kg(T, v) \ rii n c dz', 

J z 



(2) 



where i?* is the radius of the star, n; is the ion number 
density, n c is the electron number density, and T is the 
isothermal temperature. Coordinates normalized to the stel- 
lar radius are used, with (p, <f>, z) for cylindrical coordi- 
nates and (x, 9, (f>) for spherical coordinates. The factor 
K s is given by (Cox 2000): 

K s = 3.692xl0 8 (l - e - hv ' kT ^j Z^g v T-^ 2 u- 3 ,(?) 

where g v is the Gaunt factor and Z\ is the root mean 
square ion charge. For the Gaunt factor at IR wavelength, 
a power-law in frequency is adopted as given by Carciofi 
& Bjorkman (2006), 



9v = go {v/vo)~ 



(4) 



with u « 0.23 at a frequency v corresponding to 1 ^m. 
Note that the number density quantities relate to the mass 
density in terms of mean molecular weights per ion /ij 
and per free electron fj, e , respectively. Hence /^m^n; = 
p c m H n c = p. 

The next step is to consider the line transfer, for which 
the Sobolev approximation is used. Although the wind is 
in constant radial expansion, a line-of-sight gradient in 
the projected line-of-sight velocity occurs owing to the di- 
vergent expansion of the spherical wind. Since highly su- 
personic wind speeds are of interest, the Sobolev approx- 
imation remains valid except for the extreme line wings 
where spatially the velocity gradient vanishes along the 
sightline from the observer to the star center. 

For constant radial expansion, the observed Doppler 
shift in velocity units along the z-axis is 



Sobolev theory does more than simply identify isove- 
locity zones. It also offers a solution to the line transfer 
(e.g., Castor 1970; Rybicki & Hummer 1978; Hummer 
& Rybicki 1983, 1985). It is known that Sobolev the- 
ory for modeling line profiles under the assumption of 
a smooth hot star wind is inadequate for producing re- 
alistic fits to some spectral features, such as the "black 
troughs" of ultraviolet P Cygni lines (Prinja, Barlow, & 
Howarth 1990). The problem is that not only are hot star 
winds clumpy in density, but the radial profile of the wind 
velocity is non-monotonic (Lucy & Solomon 1970; Lucy 
1982). As a practical matter of line profile synthesis, the 
radiative transfer has been addressed by some in the form 
of a modified Sobolev approach, such as the SEI method 
of Lamers, Cerruti-Sola, & Perinotto (1987) or the co- 
moving frame methods of Hillier & Miller ( 1 998) or Grafener, 
Koesterke, & Hamann (2002). The non-monotonic be- 
havior is approximated as an average smooth wind with 
substantial "turbulent" broadening (e.g., Groenewegen & 
Lamers 1989). 

Since the Sobolev approximation is related to veloc- 
ity gradients, it may be inappropriate for wind clumps if 
gradients are absent throughout such structures. However, 
the situation may not be so grave for the Sobolev approx- 
imation. First, a good model for the 3-dimensional struc- 
ture (density and velocity) is currently lacking, and so the 
absence of velocity gradients remains a question. Second, 
in the scenario of this paper, the flow is treated as being 
in constant expansion, and the velocity gradient reduces 
to a geometrical one (see below) owing to spherical di- 
vergence, and clumps can certainly be expected to par- 
ticipate in the overall bulk radial expansion. Finally, the 
Sobolev approximation should remain valid for clumps 
that are small compared to the Sobolev length (i.e., L ~ 
«br/ (dv z /dr) for v^ r the dominant local broadening pro- 
cess) and that are optically thin (in contrast to treatments 
of thick clumps as in Oskinova et al. 2007). Although lim- 
itations for the line transfer in the Sobolev approach are 
noted, Sobolev is employed to the restricted problem at 
hand to explore the influence of clumping on emission 
lines and continuum formation at long wavelengths. 

The key parameter for describing the line transfer is 
the Sobolev optical depth ts for a line with opacity at 
point (p, <f>, z); the optical depth is given by 



-V~a COS 9 



(5) 



where 9 is the spherical polar angle from the z-axis as 
previously noted, and fj, = cos 9. The main value of the 
Sobolev approximation is the identification of emission 
at a fixed velocity shift in an observed line profile with a 
spatial locus of points in the flow (or "isovelocity zones"). 
In this case a fixed velocity shift transforms to a conical 
surface in the spherical wind. 



TS 



kl p A R* 
\dv z /dz\ 



(6) 



where kl is the frequency integrated line opacity with 
dimensions of area per mass times frequency, Ao is the 
central wavelength of the line transition, and the denomi- 
nator is the line-of-sight velocity gradient evaluated at the 
isovelocity zone for fixed impact parameter p. The param- 
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Red 




Fig. 1 Lower: A sightline that passes through the wind. The 
two points are intersections of the sightline with two different 
isovelocity cones shown here in cross-section. The geometry is 
for a line blend, such that at fixed A, the transition with line cen- 
ter redward ("R") of A is the rearward point, and the one that 
is blueward ("B") is nearer the observer. Above is a cartoon il- 
lustration of a line blend that might be emerge from the wind, 
with solid for the total and dotted for the individual line contri- 
butions. 

eter i?* appears because z is a normalized coordinate. In 
spherical symmetry the velocity gradient becomes 

dv z dv 
dz dx 

where x = r/R*. Equation dT) reduces to (Voo/x) (1 
/i 2 ) for a wind in constant expansion. 



cos 



V ■ 2 

- sin 

x 



(7) 



2.2 The Solution to the Radiative Transfer 

2.2.1 Solution for General Line Source Functions 

Before embarking on the steps of the derivation, more 
parameters must be introduced to describe the properties 
of the radiation. The continuum source function is repre- 
sented by Sc, which for free-free processes is Planckian, 
hence Sc = B V (T). The line source function will be Si,. 

The derivation largely follows the notation of Ignace 
et al. (2003) who elaborated on the steps originally de- 
scribed by Hillier et al. (1983). For S L = S c = B V {T), 
the total emergent intensity of continuum and line radia- 
tion is given by (eq. [Al] of Ignace et al.): 



I V =B V 1 - e -(w+rs) 



(8) 



where t$ is evaluated at the point (p, z), and r max is the 
maximum free-free optical depth along a sightline. The 
assumption of an isothermal wind is implicit in equa- 
tion © so that B v may be factored out of the integral 
to obtain a solution to the transfer equation. 



As in Ignace et al., it is convenient to express contin- 
uum and line optical depths in terms of scale parameters 
re and tl that encompass the various physical properties 
of the wind and opacity that do not vary with spatial lo- 
cation (such as T, v, etc). The Sobolev and optical depth 
can be expressed as: 



t~s = t l x 3 (1 - /i 2 ) = t l sin 8 p~ 



(9) 



which uses the fact that p = x sin 8. The continuum free- 
free optical depth from the observer to any point z in the 
wind along a fixed sightline of impact parameter p derives 
from an integral along the path as given by 

dz 



Tff = T C 



x 4 w 2 (x) 



(10) 



Tff = — r (9 — sin# cos 9) 



(ID 



Eliminating x again in favor of p and 8, the integral for the 
continuum optical depth has an analytic solution given by 

2p 3 

where tan 8 = p/z. The parameter rc is a scaling param- 
eter for the free-free optical depth that incorporates the 
various wind and atomic constants involved along with 
the frequency dependence; in relation to equation ©, it 
is given by: 



t c = R* Kg(T, v) rii.o n ,o, 



(12) 



where n^o and n c .o are the ion and electron number den- 
sities at the wind base R„. The maximum optical depth 
T max along a sightline that does not intercept the stellar 
photosphere is 



7"max — ~ g" • 



2p 3 



(13) 



When Tff 3> 1, an IR/radio pseudo-photosphere forms 
that is much larger than the hydrostatic stellar photosphere. 
In this case there is a standard integral that allows for the 
integration of the emergent intensity with impact parame- 
ter to derive the total flux of line and continuum radiation 
(Hillier et al. 1983). The new considerations include line 
blends, an additional power-law dependence to the opac- 
ity, and line source functions that need not be Planckian. 
Although a closed form for the emergent intensity can be 
given for this new generalization, there is unfortunately 
no analytic solution for the total flux, except in the spe- 
cial case of LTE for the line source functions. (That case 
will be considered in subsequent sections.) 

In deriving the emergent intensity, the first of the new 
considerations is to allow for a line blend. Suppose there 
is a wavelength where two separate lines contribute to the 
opacity. The two lines shall be referenced in terms of re- 
spective line center wavelengths as the "red" line (cen- 
tered at A/?) and the "blue" line (centered at Ab) relative 
to the wavelength A under consideration. As illustrated in 
Fig. [T] the velocity shift for the red line is necessarily at 
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an isovelocity zone that lies rearward of that for the blue 
line. Figure [TJ shows a sightline as it might cross the two 
isovelocity zones (and no more!) for the two respective 
lines at wavelength A. The following equation describes 
contributions from alternating emission and attenuations 
by the respective line and continuum opacities in the dif- 
ferent segments of the figure as working from left to right: 



h = B v (e 



Sr (1 - e 
-B v (e~ TB 



-S B U-e e 



-B v (1 - 



,-TS,B 1 



(14) 



There is an identifiable pattern that can be extended to 
allow for blends of an arbitrary number of lines. It is 
convenient to split the emergent intensity into two parts: 
the continuum contribution Ic arising from terms with 
the continuum source function and II involving contri- 
butions from all of the line source functions. 

The emergent intensity from the free-free continuum 
emission is given by 



Ic = B v 



N+l I i<j 

i - e ~ Ti ex p - rs >- 

j=l \ i=0 




(15) 



where N is the total number of lines contributing to the 
blend at a fixed frequency in the spectrum. In this for- 
mulation Tj is the free-free continuum optical depth from 
the observer to the i th resonance point along a given ray 
through the wind. In this notation tm+i = T max . Then 
Ts y i is the Sobolev optical depth at the i th resonance point, 
with rs,o = 0. 

The contribution of the total line intensity is given by 



N 




is only required when evaluating the total flux of emis- 
sion. Given that B u is constant, the above relations result 
only from a consideration of sequencing. 

The total emergent intensity along the ray is simply 
I v = Ic + It,. In order to obtain the continuum spectral 
distribution and line profile shapes, additional simplify- 
ing assumptions of constant expansion and Slj = B v 
are introduced next as a particular realization of possible 
wind models. 

2.2.2 The Solution for LTE Source Functions 

It is well-known that the level populations of atomic species 
in stellar winds is highly NLTE (Auer & Mihalas 1972). 
However, the application being considered here is for IR/radio 
spectra involving H I or He II recombination lines and free- 
free continua. The continuum that is formed at large ra- 
dius and thus roughly coincident with the line formation 
is LTE. The lines themselves are for transitions involv- 
ing high n-level values whose populations should be gov- 
erned largely by recombination and photoionization pro- 
cesses. These lines are formed at sufficiently large radius 
that in the case of the WR winds, the EUV radiation from 
the star is (a) highly diluted and (b) strongly absorbed by 
the dense wind. In this case the photoionization from the 
spatially local IR continuum that is Planckian (although 
not isotropic) will contribute to the populations of the up- 
per levels. Indeed, Hillier et al. (1983) found that the rel- 
ative upper level populations for some high n-levels of 
He II in EZ CMa were consistent with LTE. 

Here LTE is invoked as a convenience to allow for 
an analytic solution to the line and continuum fluxes and 
quick exploration of the model parameters. There is a 
tremendous simplication for the emergent continuum and 
line intensities, because all but two terms cancel exactly. 
After simplifying expression (fl4l using 5r — Sb = B v , 
the emergent intensity reduces to just 



Note that at this point, the only requirement on the solu- 
tion for equations (fl~5T > and ( TT6b is that the wind be isother- 
mal; velocity structure, ionization gradients, or variable 
clumping factors are allowed in the set of continuum op- 
tical depths Tj, The set of line optical depths rg.i are also 
quite general since they are locally evaluated only at cross- 
ings of the sightline with isovelocity zones. Indeed, it is 
not even necessary at this point to specify the velocity law 
or the geometry of the isovelocity zones; that information 



Iv — B v 



^ _ p — (Tmax + TS.R + TS,E 



B u 1 - e 



(17) 
(18) 



where tt is a sum of all three optical depths. As be- 
fore, this can be extended for an arbitrary number of line 
blends, now with tt = r max + rs,j. Each Sobolev 
optical depth must be evaluated for its respective isove- 
locity zone at a fixed frequency. This amounts to different 
velocity shifts for lines of different wavelength centers 
{A;}. These shifts can be written in normalized form as 
Wi = Vz^/voc = —fj,i. It is the total flux of emission 
that is required to simulate observables, and the intensity 
is now of a form that can be solved analytically because 
both the line and continuum optical depths have the same 
radial dependence. 
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Another important extension to previous treatments is 
that an additional power-law dependence of the opacity 
can be included. The winds are known to be clumped, and 
so this new power-law dependence can be used to repre- 
sent a radius dependence of the clumping. The clumping 
factor Dc is introduced as follows: 



Dn = 



(P) 2 



As a result, the Sobolev optical depth becomes 



T L X 

tlP 



'(3+™) (1 _ M 2 } -1 

-< 3+m ) (sine) 1+m , 



(19) 



(20) 
(21) 



where tl is a scale factor for the optical depth. Equa- 
tion ( T2lT i in terms of p and 9 will prove useful in evaluat- 
ing the total flux. The free-free optical depth becomes 



(22) 





x~ 


(4+m) dz 




— oo 










= T C p 


-(3+m) 


r 

I ( sm 






Jo 


= T C p 


-(3+m) 


x 7(m). 



\2+m 



(10 



where tc from equation ( fl2b is now expressed as 

K S {T, v)M 2 



167T 2 Rl V 2 ^ jLt e Mi TO^ 



(23) 

(24) 
(25) 



and 7 is a constant that depends solely on the power-law 
index m via the integral factor of equation (l23l . In the 
case that m = 0, 7 = tt/2. Formally, the latter expression 
only applies for rays that do not intersect the star, hence 
only forp > 1. 

As discussed by Wright & Barlow (1975), as long 
as the continuum optical depth is quite large, a pseudo- 
photosphere is formed that is signficantly larger in cross- 
section than the star itself. In this case ignoring rays with 
p < 1 represents a small error in the flux. Under these cir- 
cumstances, the solution for the total flux of emission in 
the continuum and lines is analytic. Even with line blends 
and the additional power-law factor, the form of the in- 
tegral is similar to that of Ignace et al. (2003). Here the 
steps for deriving f u are reviewed. 

The flux is given by an integral of the intesity over all 
rays: 



fu 



R 2 f 
27T ~df J ^^PdP' 



(26) 



where the integral accounts for only those spatial points in 
the wind that contribute to emission for some frequency in 
the spectrum v. Unlike the single line case considered in 
Ignace et al., this frequency point may have contributions 



from multiple line opacity sources. Using equation ( [T8l 
with constant B v , the flux becomes 



d 2 



(l-e~ TT(p) ) pdp, 



(27) 



where d is the distance to the star, and from equations (I2TI 1 
and (l24l the total optical depth is now 



TT (p) = p -(3+"0 x 

t c 7(m) + ^ T ^,i (sin6»i) 1_ 



(28) 



Again, the use of the form of the total optical depth above 
in the preceding integral relation for the flux assumes that 
the flux contribution for rays with p < 1 is relatively 
small. This is a good approximation at long wavelengths 
because the effective photosphere can indeed be much 
larger than The advantage of this treatment is that 
the integral is analytic. 

To arrive at the solution, it is useful to introduce a 



change of variable with [/( 3 +™ 1 ) — tcp ( 3+m ) 
the integral expression to the following form: 



to convert 



fu 



2/ (3+m) 



{l-cxp[£/ 3+m ( 7 

N 

-^c 1 ^TL, i sin m + 1 ^)]}, 



U^dU x 



(29) 



i=0 



where the sum refers to any of the lines contributing to 
the opacity at the frequency of interest. For N = there 
are no contributing line opacities, tl.o = 0, and expres- 
sions reduces to an expression for the continuum emission 
alone. 

The integral in that form of ( 1291 ) has a solution given 
in Gradshteyn & Ryzhik (2000) found in section 3.478, 
#2, with the final result for the observed flux /„ being 



f v = 2itB v -£-t c 



* 2/(3+m) 




N 

■E 

i=0 



TL,i 1 2\(™+ 1 )/ 2 
— (1—0 



2/(3+m) 



(30) 



where T is the Gamma-function, and sin 6 = (1 — w 2 ) 1 / 2 . 
For reference, the solution equation (f30b shall be referred 
to as the "Free-free + Recombination Emission line in 
Constant velocity" model, or the "FREC" model. 
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Fig. 2 Illustration of the analytic solution for a continuum dominated by free-free emission and He II recombination lines, such 
as occurs for WR stars. Spectra are shown as normalized to the continuum fc. The two plots left and right are for two different 
wind terminal speeds as indicated. The upper and lower panels are for two different power law indices m: upper is for constant 
clumping factor and lower is for clumping that decreases inversely with radius. Note that the analytic results account fully for line 
blending. The different curves within each panel are for different line strengths. 



In regions of the spectrum that are absent of lines, the 
FREC model provides the continuum flux. In the Rayleigh 
Jeans limit, the continuum reduces exactly to a power-law 
form. Assuming hv <C kT and using the Gaunt factor of 
equation (0|, the continuum distribution between lines be- 
comes a function of just one free parameter m (i.e., for the 
value of u known): 



f v OC j,( 2 + 2 '»)/(3+m) g 2/(3+m) 



2(m+l— u)/(3+m) 



Note that with u appropriate for the radio band, and m — 
0, the 0.6 power-law slope of Wright & Barlow (1975) is 
recovered. 

Returning to the lines, why does the line shape vary 
with m-value? Although both line and continuum opac- 
ities scale with the square of the density, the line optical 
depth has an additional dependence with angular escape 
through a velocity gradient. As a result, the line shape is 
not invariant with respect to the continuum formation as 
m is varied. 

Examples of model spectra produced by the FREC 
solution are shown in Fig. [2] The models are continuum- 
normalized to focus on the line shape effects and the in- 
clusion of blending. The left and right plots are for dif- 
ferent wind speeds of 2500 and 5000 km s _1 . All panels 
show a segment of the 13-14 micron spectrum. All emis- 
sion lines are of He II recombination; the strong line at 
13.12 microns is He II 10a. The upper panel is for m = 0, 
corresponding to a wind that has a constant clumping fac- 



tor. The lower panel is for m = 1, indicating a clumping 
factor that decreases with radius. 

To clarify better the line-shape dependence on the in- 
dex m, it is useful to consider limiting cases of a single 
strong or weak line: a strong line occurs for 17, 3> rc and 
a weak one for tl <C tq . A single strong line yields the 
profile shape given by 

t \-\ „„21 (l+™)/(3+m) 

(31) J» x L 1 ~ w \ ■ ^ 2 > 

In the case of a weak line, the flux is of the form f v oc 
1 + e, where e has the functional dependence of 

21 (l+m)/2 



eoc [l-v?\^"'"". (33) 
The profile shapes for thin and thick lines are not the 
same. Thick lines tend to be more "rounded" or "bubbled" 
as compared to thin ones. The half width at half maximum 
(HWHM) of a thick line is greater than for a thin one. In 
the case of a thick line, it is 



HWHM thick = Voo W 1 - 
whereas for a thin line, it is 



(3+m)/(l+m) 



HWHM thil 



2/(l+m) 



(34) 



(35) 



For m = 0, HWHM thick = 0.94 Voo and HWHM thin = 
0.87woo. The difference in half-width from the full ter- 
minal speed is small, but it becomes greater for increas- 
ing values of m. With m = 1, the half-width values be- 
come 0.87woo versus 0.71woo. For fast winds of WR stars 
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Fig. 3 Comparisons between thick and thin emission lines 
from the FREC model. The two top panels display line pro- 
files for thick (left) and thin (right) lines. These are continuum 
subtracted and then normalized to their peak values and plotted 
against normalized velocity shift. Each curve is for a different 
m value, from m = —0.5 to 4 in intervals of 0.5. At bottom is 
a plot of the HWHM of thick and thin lines against m. The line 
widths are normalized to the wind terminal speed. 



at around 2000 km s , the difference amounts to over 
300 km s -1 , which is potentially measurable. Figure [3] 
compares, as a function of m value, thick versus thin line 
profile shapes (top two panels) with a corresponding plot 
(bottom panel) of the HWHM values relative to the wind 
terminal speed. 

It is also useful to consider trends in the line equiva- 
lent width (EW). The line equivalent width for the FREC 
model with a single line is given by: 



EW 



A v c 



TL 



7 (to) tc 



(36) 
dw, 



where w = i> bs/voo ^ s trie normalized velocity shift across 
the line. Figure |4] shows curves for the EW normalized 
to XqVoc/c and plotted as a function of the relative line 
strength tl/7(to)tc- The different curves are for differ- 
ent values of m at increments of 0.25. Curves are closely 
clustered with shallow slopes for larger values of m rep- 
resenting steeper declines in the clumping factor with ra- 
dius in the wind. 



2 4 6 

Relative Line Strength 

Fig. 4 Normalized line equivalent widths (EW) are plotted 
against the relative line strength (see text). Each curve is for a 
different power-law index m, ranging from m = —0.5 up to 
m = 3 in increments of 0.25. The EW is seen to grow in line 
strength as expected, but the rapidity of growth is much less for 
clumping factors that are more steeply decreasing with radius. 



The overall result is that increasing values of m steepen 
the decline in clumping with radius, narrow the line width, 
increase the relative line strength (because 7 decreases 
with m), and decrease the line EW. Equations (f2TT > and 
(l24l i show how the line and continuum optical depths change 
with to, and the rate of decline for these is not the same 
along a fixed sightline. The line becomes thinner more 
rapidly at a given resonance point as compared to the to- 
tal line-of-sight free-free optical depth. 

3 Application to WR 90 

With the solution to the radiative transfer through the wind 
at hand, it is useful to consider an illustrative model. The 
continuum slope is set by equation OTb . To construct a 
full emission line spectrum, a tabulation of line optical 
depth scales {r^ i} are required as input to the model. The 
IR spectra of WR stars are characteristically He II emis- 
sion line spectra, and at wavelengths longer than about 10 
microns (relevant to the Spitzer/TRS), the line transitions 
involve somewhat high n-values, levels above n = 9. 

With the LTE assumption, the Sobolev optical depth 
scale factors vary with the line transition as follows: 



t l (xgfX 



1L 

fJ.i 



(37) 
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where / is the oscillator strength, g values are statisti- 
cal weights for respective levels, and n is a scaled num- 
ber density (i.e., no dependence on radius). Assuming the 
Rayleigh-Jeans limit as a simple case, the final result of 
this is that 

tl oc gf oc A 2 Aji n 3 -. (38) 
9j 

It is now easy to compare the line optical scale factors for 
any two lines tx and t' l involving different upper levels j 
and j'. Given in ratio, the optical depth scales are related 
as follows, 




A FREC spectral model depends on just three items: 
(1) m that sets the continuum slope, (2) To a fiducial line 
optical depth for some line from which all other lines 
scale, and (3) a list of He II line data. An example appli- 
cation in Fig. [5] shows a FREC model overplotted with 
Spitzer IRS/SH data of the WC star, WR 90. Data are 
shown as solid, and the model is the dotted curve. Tran- 
sitions for Hell up to n = 30 are included. The ob- 
served spectrum includes lines of other atomic species be- 
sides just recombination lines of He II, notably the boom- 
ing line of [Nelll] at 15.55 microns, but also lines from 
ionized carbon. The model makes no attempt to include 
any lines other than He II recombination lines. The dot- 
ted curve is the simple two-parameter spectral model as 
"eyeballed" for a good fit. The model does adopt the wind 
terminal speed of WR 90 from Ignace et al. (2007) and 
has been gaussian convolved to match the resolution of 
the IRS in SH mode. 

A value of m — 0.9 reasonably reproduces the ob- 
served slope. After fixing the continuum slope, the optical 
depth scale was varied until approximate matches to the 
observed emission in Hell 10a at 13.12 /im and 11a at 
17.26 fim were achieved. The fit is only illustrative, and 
no claim is made that m = 0.9 is an accurate descrip- 
tion of the clumping distribution in the wind of WR 90. 
Modeling has shown that the winds of WC stars have 
rather extended acceleration zones to large radii (Hillier 
& Miller 1999; Grafener & Hamann 2005), hence the un- 
derlying assumption of constant radial expansion in the 
region where the 10-20 micron emission forms may not 
be valid for WR 90. Still, this example application to 
Spitzer data shows that the FREC model can give qual- 
itatively reasonable fits. 

4 Summary 

In this paper an elegant result originally derived by Hillier 
et al. (1983) has been expanded significantly to allow for 
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Fig. 5 Shown is Spitzer data for the WC star, WR 90 (solid), 
along with a model fit using the FREC model (dotted). This 
is a two-parameter fit: m = 0.9 is chosen to match the con- 
tinuum slope, and the lines are governed by a single optical 
depth parameter assuming LTE for relative level populations. 
The spectrum has been gaussian convolved to the resolution of 
the IRS/SH instrument. The stronger lines of He II are approxi- 
mately matched. Note the prominence of [Ne III] 15.55 microns, 
which is not part of the FREC and would need to be modeled 
separately (e.g., Smith & Houck 2005). 

greater versatility. Although the assumptions remain re- 
strictive (spherical symmetry, isothermal, and the con- 
stant radial expansion), the solution to the radiative trans- 
fer is analytic and allows for both line blends and a power- 
law dependence of the clumping with radius. The latter is 
especially significant since clumping is now understood 
to be an important consideration in deriving mass-loss 
rates from observations and since there are outstanding 
questions regarding the evolution of the clumping factor 
throughout stellar winds. 

What is especially interesting about the new FREC 
models is a relation between line widths and line equiv- 
alent widths that has a dependence on the exponent m. 
Quantitatively the predicted relations may not be realized 
in data owing to the various assumptions built into mak- 
ing the FREC solution analytic; however, models based 
on the generalized line source functions from equation ([Tol l 
may be new avenues of consideration for quick parame- 
ter exploration for those who use the more sophisticated 
radiative transfer techniques for spectral synthesis of stel- 
lar winds. Archival data from the Infrared Space Obser- 
vatory or the Spitzer satellite may have adequate spec- 
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tral resolution for sources with fast winds to test against 
such relations and place limits on the variation of clump- 
ing with radius. 

Although this paper focused on Hell emission line 
spectra of WR winds for application, owing to the high 
densities of these winds, the solution may be applicable 
to other sources, like the hydrogen recombination line 
spectra of Luminous Blue Variable (LBV) stars, such as 
P Cygni (Lamers et al. 1996). The FREC solution may 
even have application data to suitable hypercompact or ul- 
tracompact H II regions. These sources show radio recom- 
bination lines and free-free continuum spectra plus evi- 
dence for density gradients (e.g., Keto, Qizhou, & Kurtz 
2008). 
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